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Simulations of vortex tube dynamics reveal that the non-Gaussian nature of 
turbulent fluctuation originates in the effect of random advection. A similar non- 
Gaussian distribution is found numerically in a simplified statistical model of random 
advection. An analytical solution is obtained in the mean-field case. 
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In the study of statistical physics for non-equilibrium systems the deviation from the 
Gaussian distribution has been the central issue. We are expecting the existence of universal 
mechanisms of producing non-Gaussian distributions as in the case of thermal equilibrium 
systems, however, our knowledge has been still in the elementary level. 

Non-Gaussian probability distribution function (NGPDF) is especially important in fluid 
turbulence. Turbulence is a typical far-from-equilibrium phenomena as it is characterized 
by the energy cascade which violates both the detailed balance and equi-partition. The 
appearance of NGPDF is so common and we can not construct any theory without taking the 
non-Gaussian nature into account |TJ.0]. Latest technique of direct observation of turbulence 
H and also direct numerical integration of Navier-Stokes equation || are clarifying the 
details NGPDF from experimental viewpoints. 

Recently, theoretical investigation to NGPDF in turbulence is attracting much attention. 
An approach is to assume the multifractality in the geometry of turbulent velocity field as 
an ideal limit |Q, and another one called the mapping closure conjectures the existence of a 
kind of smooth map which describes the time evolution of statistical quantities |5|||. 

Yakhot et al. JF| also present a phenomenological excellent approach in order to derive 
a NGPDF for vorticities in the fluid turbulence. Although it is fairly reasonable, physical 
meanings in their assumptions are not clear. 

In this paper we first introduce a numerical model of vortex dynamics on lattice and show 
that the vorticity distribution in a turbulent steady state is actually far from Gaussian but 
closer to a symmetric exponential distribution. In order to clarify the origin of this NGPDF 
we modify the dynamics so that vortex tubes move randomly by passive advection. It is 
shown that this modification does not seem to affect the vorticity distribution, which implies 
that the random advection is playing an important role. Then, we focus our attention to 
the effect of random advection and introduce a random scalar advection model on lattice. 
Surprisingly the scalar also follows a distribution which is almost identical to the former 
NGPDFs. A mean-field version of the scalar model is then solved analytically showing a 
clear evidence of NGPDF. 

Let us first introduce our lattice vortex model || which is equivalent to the following set 
of vorticity equations of incompressible fluid in the continuum limit 

(uj ■ V)v + z/AcD, (la) 

uj, (lb) 
0, (lc) 



duo 
~dt 



+ (v-V)uj 

YOtV 

divv 



where v, uj and v are the velocity, vorticity and viscosity, respectively. We assign vor- 
ticities on a simple cubic lattice. Every bond has only one vorticity component along its 
direction. This means that an x-bond has only x component of vorticity, i.e., a bond can 
be viewed as a vortex tube. For a given configuration of vorticities the vorticity field is 
calculated by using the Biot-Savier law (equivalent to eqs.(|T5D and (|lc|). The fluid velocity 
of a bond at r,v(r), is estimated at the mid-point of the bond. Dynamics is defined so as 
to satisfy the vorticity equation( 



Due to the advection term (v ■ V)u; in eq. (|TaD the vorticity flows with the fluid as shown 
by the name of Kelvin's theorem ||10|| , the vorticity uj z at a z-bond is transported to its six 
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neighboring z-bonds. By this effect the value of uj z at r ± i is increased by ±Jj 2 At/2, where 
i is a unit vector directing either x, y or z direction, and Jy is the flux of vorticity, ViUJj. 
The coefficient 1/2 of the vorticity flux and the signs ± are introduced to keep the spatial 
symmetry. The first term on the right hand side of eq.(|Ta|) shows the vortex stretching 
term, so we have to modify the vorticities around the advectively transported bonds so that 
modified vorticities make loops as shown in Fig.|I]. In this procedure we add ^fJi Z At/2 to the 
bonds at r ± i/2 ± z/2, namely, x and y component appear due to the elongation of vortex 
tube. The diffusion term in eq. (|Ta| ) is discretized by the usual finite difference method. 
It has been shown that this dynamics is equivalent to eq.flla]) in the continuum limit and 
numerical simulation of this vortex tube dynamics on a periodic 24 x 24 x 24 lattice has 
reproduced many basic properties of fluid turbulence ||. 

As for the probability distributions we obtained the following results: 

1. The velocity components (for example v x ) follows nearly Gaussian. 

2. The distribution of differentiated quantities, such as dv x /dx,dv x /dy and u> x , are 
non-Gaussian and close to exponential (see FigfJ). 

3. By removing lower wave number components, the distribution of v x also becomes 
closer to an exponential as firstly found by She et al JTT|. 



We also observe a kind of conditional distributions for velocities on bonds whose vortici- 
ties are in a fixed range. The distribution are nearly identical both for large u and for small 
u, indicating that the velocity on a bond is nearly independent of its vorticity strength. 

We now modify our model to seek the origin of the NGPDF. The modification is to 
randomize the velocity field at each time step keeping its distribution. Namely, we do not 
use Biot-Savier law, but the vorticities are transported and elongated by random advec- 
tion. The probability distribution of uo after several hundreds time steps is quite similar to 
that of original lattice vortex model (see Fig||). This result clearly shows that the appear- 
ance of NGPDF is independent of the details of velocity field. As suggested by Sinai and 
Yakhot |14]J7[] NGPDF may be caused by random advection, and if so, we can expect a large 
universality class including another exponential-like distribution in thermal convection flow 

0,0 ■ 

In order to see the effect of random advection more clearly we introduce a scalar model 
on lattice. The model is governed by the following stochastic equation 

w(r , t + At)= u(r , t) - v(r Q , t)co(r , t) At + £ P(r , r, t)v(r, t)u(r, t)At, (2) 

r 

where v(r, t) and P(r , r, t) are independent random variables. v(r, t) takes non-negative 
values, and P(r , r, t) — 1 when a flow from site r to r exists at time t, otherwise P(ro, r, t) = 
and it is normalized as J2 r P{ r o, r , t) = 1. The one- dimensional version of eq.(@) is defined 
by the special case that either P(r , r — 1, t) or P(r , r + 1, t) is equal to 1 with probability 
1/2. We perform the simulation on a 1-dimensional lattice of size 10000 with the periodic 
boundary condition. The random number v (r, t) is in the range of [0, 0.5] and At = 0.5. For 
a given random initial condition of u(r, 0) we observe the distribution of u(r, t) at sufficiently 
large t. The result is also plotted in Fig.|2|. The distribution is far from Gaussian and is very 
close to those of former cases. Now it seems obvious that the origin of NGPDF is in the 
simple random advection transport. 

A mean-field version of eq. (0) can be solved analytically. We consider the situation that 
a site is interacting with a mean-field site. To avoid unimportant complexity we assume the 
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case that the random number v takes either or jo/ At with probability 1/2 (jo G [0, 1]). In 
this situation the stochastic eq.(Q) becomes 

u(t), Prob.1/4 

i u{t)+3 u M , Prob.1/4' v ; 

(1 - j )u(t) + joU M , Prob.1/4 

where ujm is an independent random number having the same distribution as u(t). By 
introducing the characteristic function 



/oo 
e ipuJ P(cu,t)duj, 
-oo 



eq-(l) b 



ecomes 



Y(p, t + At) = ^{Y(p, t) + Y(p - pjo, t)}{l + Y(p 3o , t)}. (4) 

Taylor expansion of eq. (§) in terms of p gives a set of equations for the moment functions 
{< u n >}, and it can be easily shown that a steady state exists when < uj >^ 0. 

The steady state solution is obtained with the aid of algebraic calculation by computer. In 
Fig.|3] we plot the cumulants for orders up to 7 together with those for one-sided exponential 
distribution for comparison. Higher order cumulants are not zero in any case and showing a 
tendency to diverge as the order goes to infinity, which clearly demonstrates that the steady 
state distributions are not Gaussian. 

In the special case of jo = 1/2 we can estimate the functional form of Y(p) as follows. 
By denoting Y n = Y(p2 n ) eq.(§J) becomes 

Y n +i — Y n = 2 • ™ — — . (5) 

Supposing that n can be a continuous number and we approximate the left hand side by 
a derivative dY/dn. Then eq.([D can easily be integrates and we have the following cubic 
equation for Y(p) with < u >— 1, 

p 2 Y(pf + (Y(p) - l) 2 = 0, (6) 

which can be solved explicitly by Cardano's formula. The asymptotic behavior for p — > oo 
is given as 

Y(p) cc p-^. (7) 

This power law decay is very different from the case of a Gaussian distribution where 
the characteristic function decays faster than any power of p. In the case of exponential 
distribution the characteristic function is given by 1/(1 + ip) (asymmetric case) or 1/(1 + 
p 2 ) (symmetric case), which also shows a power law decay as p —>■ oo. In this sense our 
distribution is far from Gaussian, but has a similarity to the exponential distribution. 

For smaller j the cumulants are smaller, and in the limit of jo — > it can be shown that 
the normalized cumulant of order n, < cu n > c / < to 2 >^ 2 , vanishes for n > 3. Namely, the 
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distribution converges to a Gaussian in this limit, which agrees with the well-known fact 
that the distribution function satisfying the usual diffusion equation in the continuum limit 
is a Gaussian since eq.(f^) becomes a usual diffusion equation in the limit of At — > 0. This 
indicates that the finiteness of At is very important for the appearance of NGPDF in our 
random advection dynamics. 

One may think that these results for At 7^ are meaningless because we must take 
the limit of At — > in order to recover eq. (|Ta|) . This is a wrong argument neglecting the 
finiteness of correlation time in the real dynamics with At = 0. It is known that the fastest 
mode in real turbulent flows has a characteristic time approximately given by yu/e, where 



e is the mean energy dissipation rate [jnj]. For time scales less than this value advections 
can not be treated as random. In our random advection models At should be viewed as the 
correlation time in real dynamics, therefore, it should take a finite value. 

Our theoretical approach is successful only in the case of < uj >^ 0, while observed 
distributions are symmetric with < uj >= 0. At present we have no rigorous way of con- 
necting these two cases, but we are now considering the effect of spatial fluctuation as the 
key of solving this difficulty. Even in the case of < uj >= averaged over the whole space, 
< uj > can be nonzero if the average is taken over a finite area due to the spatial fluctuation. 
So, if the distribution of uj is determined rather locally not using the information from the 
whole space, then the distribution for < uj >^ may have direct consequence to the real 
distributions. 

In summary, we showed that the random advection dynamics creates far-from Gaussian 
fluctuations whose distribution functions are closer to the exponential. As shown in the dis- 
cussions of the statistical scalar model on lattice and its mean-field analysis, the mechanism 
of NGPDF is so simple that a wide application may be expected. 
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FIGURES 



FIG. 1. Dynamics of vortex tubes. Direction of arrows indicates the sign. (See text) 

FIG. 2. Probability distributions of ui x ; O for true vortex dynamics H, + for random advection 
model of vortex tubes, □ for the 1-dimensional scalar model. They are averaged over ten samples. 
Error bars stand for standard deviations. uj x is normalized so as to have variance of unity. 

FIG. 3. Higher order cumulants {oo n ) c . + for the mean-field model; jo = 0.1,0.7 and 0.9 from 
bottom to top. The exponential distribution is shown for comparison(O). 
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